clc;
A=[1,2,3;
   4,5,6;
   7,8,0];
C=A^(1/3),
e=norm(A-C^3)


j1=exp(sqrt(-1)*2*pi/3);

A1=C*j1,A2=C*j1^2

e1=norm(A1-A1^3),e2=norm(A-A2^3)%norm用来计算矩阵范数

A=sym([1,2,3;4,5,6;7,8,0]);
C=A^(sym(1/3));
C=vpa(C);
norm(C^3-A)%实际上这个地方是无法运行成功的
